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Abstract 

Classical Pitowsky correlation polytopes are reviewed with particular emphasis 
on the Minkowski- Weyl representation theorem. The inequalities representing the 
faces of polytopes are Boole's "conditions of possible experience." Many of these 
inequalities have been discussed in the context of Bell's inequalities. We introduce 
CddlF, a Mathematica package created as an interface between Mathematica and 
the ccld program by Komei Fukuda, which represents a highly efficient method to 
solve the hull problem for general classical correlation polytopes. 
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1 Boole-Bell Type Inequalities And Their Geometric 
Representation 

In the middle of the 19th century the English mathematician George Boole formulated 
a theory of "conditions of possible experience" [[j], ||, ||, ^, ||]. These conditions are re- 
lated to relative frequencies of logically connected events and are expressed by certain 
equations or inequalities. More recently, similar equations for a particular setup which 
are relevant in the quantum mechanical context have been discussed by Clauser and 
Home and others ^ [A |h, Pitowsky has given a geometrical interpretation in terms of 
correlation poly topes [K], [|. 

1.1 Simple urn model 

Consider an urn containing some balls of different colors and styles: Each ball can be 
described by two properties, let us say "yellow" and "wooden", so each ball can have 
the property "yellow" or the property "wooden", but it can also have both "yellow and 
wooden". The state of the urn can be given by the probabilities to draw a ball with one 
of these properties: pi is the proportion of yellow balls in the urn, pi the proportion of 
wooden ones and pn denotes the proportion of yellow and wooden balls. If there are 
enough balls in the urn these proportions are in fact the probabilities to get a ball with 
the special property by drawing. Clearly the inequalities 

0<pn< P 2<\ and 0<pn<p\<l (1) 

are fulfilled by the proportions and so pi, p2 and pn can be seen as probabilities of 
two events and their joint event only if these inequalities are satisfied. Simply by taking 
some appropriate numbers (pi = 0.6, p2=0.72 and pn=032) we can see, that equations 
(|l]) are not sufficient. If we take the probability to draw a ball which is either yellow 
or wooden (pi + p% - pn) into consideration, a new inequality can be found that is not 
satisfied by the numbers chosen: 

0<pi+p 2 -pi2< 1 (2) 

It can be shown that the inequalities ([[]) and (fy are necessary and sufficient for the 
numbers p\, pi and pn to represent probabilities of two events and their joint 

1.2 Geometrical interpretation 

Itamar Pitowsky g f§ |(||] has suggested a geometric interpretation. Consider the 
truth table [j]of the above urn model, in which ai and represent the statements that 
"the ball drawn from the urn is yellow," "the ball drawn from the urn is wooden," and 
in which an represent the statement that "the ball drawn from the urn is yellow and 
wooden." The third "component bit" of the vector is a function of the first components. 
Actually, the function is a multiplication, since we are dealing with the classical logical 
"and" operation here. Let us take the set of all numbers (pi, pi, pn) satisfying the 
inequalities stated above as a set of vectors in a three-dimensional real space. This 
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Table 1: Truth table for two propositions ai,a? and their joint proposition a 12 = a\ A«2 



amounts to interpreting the rows of the truth table as vectors; the entries of the rows 
being the vector components. This procedure yields a closed convex poly tope with 
vertices (0,0,0), (1,0,0), (0,1,0) and (1,1,1) (cf. Figure [l]). The extreme points (vertices) 
can be interpreted as follows: 

(0,0,0) is a case where no yellow and no wooden balls are in the urn at all, 
(1,0,0) is representing the configuration that all balls are yellow and no one is wooden. 
(0,1,0) is representing the configuration that all balls are wooden and no one is yellow. 
(1,1,1) is a case with only yellow and at the same time wooden balls. 




Figure 1 : Polytope associtated with the urn model 



1.3 Minkowski- Weyl representation theorem 

The Minkowski-Weyl representation theorem (e.g., flilj , p. 29]) states that compact 
convex sets are "spanned" by their extreme points; and furthermore that the represen- 
tation of this polytope by the inequalities corresponding to the planes of their faces is 
an equivalent one. 

Stated differently, every convex polytope in an Euclidean space has a dual descrip- 
tion: either as the convex hull of its vertices (V-representation), or as the intersection of 
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a finite number of half-spaces, each one given by a linear inequality (H-representation) 
This equivalence is known as the Weyl-Minkowski theorem. 

The problem to obtain all inequalities from the vertices of a convex polytope is 
known as the hull problem. One solution strategy is the Double Description Method 
JT^ ] which we shall use but not review here. 



1.4 From vertices to inequalities 

For the above simple urn model, the inequalities are rather intuitive and can be easily 
obtained by guessing. This is impossible in the general case involving more events and 
more joint probabilities thereof. In order to obtain the relevant inequalities — Boole's 
"conditions of possible experience" — we have to find a hopefully constructive way to 
derive them. 

Recall that a vector is an element of the polytope if and only if it can be represented 
as a certain bounded convex combination, i.e., a bounded linear span, of the vertices. 
More precisely, let us denote the convex hull con\{K) of a finite set of points K — 
{x 1 ,...,x n }e^by 

conv(7r) = |A,iXi+ - +A, n x n n> > 0,]^ = 1 j . (3) 

In the probabilistic context, the coefficients X, are interpreted as the probability that the 
event represented by the extreme point xj occurs, whereby K represents the complete 
set of all atoms of a Boolean algebra. The geometric interpretation of K is the set of all 
extreme points of the correlation polytope. 

In summary, the connection between the convex hull of the extreme points of a 
correlation polytope and the inequalities representing its faces is guaranteed by the 
Minkowski-Weyl representation theorem. A constructive solution of the corresponding 
hull problem exists (but is NP-hard [0). 

For the special urn model introduced above this means that any three numbers (p\, 
P2 and P12) must fulfill an equation dictated by Kolmogorov's probability axioms 

(pi,P2,Pn) = U (0,0, 0)+X2(0,l,0)+A-3(l,0 ) 0)+X 4 (l, 1,1) = Qv+U,)* + U,U)- 

(4) 

It is important to realize that these logical possibilities are exhaustive. By definition, 
there cannot be any other classical case which is not already included in the above 
possibilities (0,0,0), (1,0,0), (0, 1,0), (1, 1, 1). Indeed, if one or more cases would be 
omitted, the corresponding polytope would not be optimal; i.e., it would be embed- 
ded in the optimal one. Therefore, any "state" of a physical system represented by a 
probability distribution has to satisfy the constraint 

\l + Xz + h3 + U = l- (5) 

The four extreme cases A., = 1, Xj = for i E {1,2,3,4} and j ^ i just correspond to 
the vertices spanning the classical correlation polytope as the convex sum 

A generalization to arbitrary configurations is straightforward. To solve the hull 
problem for more general cases, an efficient algorithm has to be used. There are some 
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algorithms to solve this problem, but they run in exponential time in the number of 
events, thus it can be solved only for small enough cases to get a solution in conceivable 
time. 



1.5 From inequalities to vertices 

Conversely, a vector is an element of the convex polytope if and only if its coordinates 
satisfy a set of linear inequalities which represent the supporting hyper-planes of that 
polytope. The problem to find the extreme points (vertices) of the polytope from the 
set of inequalities is dual to the hull problem considered above. 

1.6 Quantum mechanical context 

In the quantum mechanical case the elementary irreducible events are clicks in particle 
detectors and the probabilities have to be calculated using the formalism of quantum 
mechanics. It is by no means trivial that these probabilities satisfy Eq. (Q), in particular 
if one realizes that quantum Hilbert lattices are nonboolean and have an infinite number 
of atoms. As it turns out, Boole's "conditions of possible experience" are violated if 
one considers probabilities associated with complementary events, thereby assuming 
counterfactuality. (This is a development and a generalization Boole could have hardly 
forseen!) 





Figure 2: Experimental setting to test the violation of Boole - Bell type inequalities 



As an example we take a source that produces pairs of spin-^ particles in a singlet- 
state (|\|/) = jj,) — | It)))- The particles fly apart along the z axis and after the 
particles have separated, measurements on spin components along one out of two di- 
rections are made. If, for simplicity, the measurements are made in the x-y plane 
perpendicular to the trajectory of the particles, the direction of the measurement can 
be given by angles measured from the vertical x axis (oci and 0C2 on the one side, Pi 
and P2 on the other side). On each side the measurement angle is chosen randomly 
for each pair of incoming particles and each measurement can yield two results - in | 
units: for spin up and for spin down (cf. Figure ||). 

Deploying this configuration we get probabilities to find a particle measured along 
the axis specified by the angles OCi, OC2, pi and P2 either in spin up or in spin down state 
denoted as p a \, p a 2, Pbi, Pb2 - and we also take the joint event of finding a particle on 
one side at the angle OCi (0C2) in a specific spin state and the other particle on the other 
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side along the vector pi (P2) in a specific spin state, denoted as p a ibi, Palbl, Palbl and 
Palbl- To construct the convex polytope to this experiment we build up a truth table 
of all possible events using a "1" as "spin up is detected along the specific axis" and a 
"0" as "spin down is detected along the specific axis" (table ||). The rows of this table 
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Table 2: Truth table for four propositions 



are then identified with the vertices of the convex polytope. By using the Minkowski- 
Weyl theorem and by solving the hull problem, the vertices determine the hyper-planes 
confining the polytope, i.e. the inequalities which the probabilities have to satisfy. As 
a result the following inequalities are gained: 



< Paibi < Pai < 

b,<l 



Pai + Pbi — Paibi < 



i<l,0<Paibi<Pbi<l t = l,2 

i = 1,2 



(6) 



(7) 



— 1 < Pa\b\ + Palbl + Palbl — Palbl ~ Pal ~ Pb2 < 
- 1 < Palb 1 + Palbl + Palbl - Palbl ~ Pal - Phi < 
- 1 < Palbl + Palbl + Palbl ~ Palbl ~ Pal ~ Pbi < 
- 1 < Palbl + Palbl + Palbl ~ Palbl ~ Pal ~ Pbi < 

The last four inequalities are known as Clauser-Horne inequalities. As noticed 
above the probabilities have to be seen in a quantum mechanical context. If we consider 
the singlet state of spin-^ particles \\\f) = 4=(| fj.) — | J,|)) it is well known that the 
probability to find the particles both either in spin up or in spin down states is given 
by pTT (Q) — pJ-J-(0) = jsin 2 (8/2) - where 8 is the angle between the measurement 



directions. The single event probability is clearly p, = j. By choosing 
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as measurement directions, we get for p = (p a i , p a2 , pbi , Pbi , Pa l b l , Pdib uPaWliPdaa)'- 
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and one of the inequalities found in (0) is violated: 



3 3 3 111 

Pa\b\ + Palbl + Pa2b2 ~ Palbl ~P al ~P b2= $ + % + g~°~2~2 = 8 >0 ^ 

The generalization is straightforward. Violations of certain inequalities involving 
classical probabilisties — Boole's "conditions of possible experience" [Q] — also appear 
in higher dimensions in configurations containing more particles and/or more measure- 
ment directions. We shall consider more examples below. 



2 Installation 

2.1 Mathematica 

All functions described in the following section can be found in the Mathematica- 
package cddif.m. In general this package has to be loaded into the current Mathematica- 
kernel by the command <<'path to cddif.m' /cddif.m, short description and 
usage of the functions is available by entering ? ' <f unction>' . 

To guarantee a proper run of all functions it is necessary (and hopefully suf- 
ficient) that cdd is located in any directory listed in the PATH-variable (usu- 
ally /bin, /usr/bin, /usr/local/bin, ...)[] or in the current working directory, 
which can be shown by evaluating Directory [ ] or changed using the function 
SetDirectory [ directory. String] . If one likes to avoid the frequent use of this func- 
tion one can append this command to the package-file cddif.m before the line End [ ] 
so that on each loading of the package the directory is set automatically to a personal 
working directory. 



2.2 cdd 

cdd is a C++ (ANSI C) implementation of the Double Description Method by 
Komei Fukuda[|l4]|. It generates all vertices (i.e. extreme points) and extreme rays of 
a general convex polyhedron given by a system of linear inequalities. Conversely, it 
solves the hull problem by generating a system of linear inequalities given all vertices. 

At this point we refer to the documentation of the program for the installation of the 
cdd - package, in particu lar to the file cdd. readme included in the package and to the 



online documentation at [http://www.ifor.math.ethz.ch/~fukuda/cddman/cddman.html . 



cdd is available for free and one can download it from the homepage of Komei 
Fukudajl"^ (http://www.ifor.math.ethz.ch/Jukuda/cddJiome/cdd.html), where one can 
also find a manual to the usage of cdd, especially descriptions to the format of the 
input- and output-files and to options that can be passed to cdd. 

1 for setting environment variables look at the manual of set or env 
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2.3 Installation on Windows-platforms 



Currently a version of cdd executable on Window s platforms can be down- 
loaded from frttp://tph.tuwien.ac.at/~svozil/cdd/cdd.exe (a different c ompilation is at 
[http://www.wis.kuleuven.ac.be/wis/algebra/kathleen/files/cdd06 1 .exe ). 
On Windows-systems Mathematica must be able to find cdd in a directory listed in 
the PATH-variable or in the current working directory, too. To set the PATH-variable 
in Windows 2000/NT go to the "control panel" and click on the "system properties", 
then click "advanced" and there is a place where the variable PATH is specified. Here 
one can add the path to cdd (separated by a semicolon). In WindowsME one needs 
to go execute "msconfig" to get to the System Configuration Utility - in "Environ- 
ment" one can edit the PATH-variable and for Windows98/95 one must edit the file 
"autoexec.bat" to get the path set. 

Finally it can be necessary to rename the cdd - executable file (e.g. from cdd061.exe) 
to cdd.exe. The CddlF-Package uses cdd as default command to run cdd, using 
the function SetCddCmd [ cmd_String] one can change this behavior. Like already 
stated above one can also add this commandline to the package-file cddif.m just before 
the End [ ] -statement to change the default command automatically when loading the 
package. 



3 Description Of Functions 

In this section all functions of the CddlF - package are listed. For each function the 
syntax including the necessary parameters (if parameters are optional, it has an " opt: " 
as prefix), a description and an example is given. 

3.1 CddFormat 

CddFormat [ vertices, opt: options] 

vertices (List): List of m vertices in n dimensions of the form 

{{•^1 1 r %12r • • • r X\n \ r {^21 r • • • r ^2n \ r • • • r {-^ml r Xm2 r • • • r %mn\} 

options (List): Options to cdd (e. g. adjacency, nondegen- 

erate, minindex,...) see documentat ion to cdd 

( http://www.ifor.math.ethz.ch^fukuda/cddman/cddman.html ) 

Description: A list of vertices, which can be determined for example by 
TruthTable [..., Include Vars^False] , are converted to a format recognized by cdd. 
Additionally options to cdd can be declared. 

Example: 

m[]= CddFormat [ {{1,0,0}, {0,1,0}, {1,1,1}}] 



9 



out[]= {v-representation begin, {3, 4, integer}, 
{1, 1, 0, 0}, {1, 0, 1, 0}, {1,1, 1, 1}, 
end} 

3.2 ToCddExtFile 

ToCddExtFile [ file, vertices, opt: options] or 
ToCddExtFile [ file, particles, measurements, opt: options] 

file (String): Filename for output of H-representation ( ".ext "-suffix is automatically ap- 
pended) 

vertices (List): List of m vertices in n dimensions of the form 

{{•^1 1 i ^-12 r • • • r X\n\ , {"£21 r • • • r %2n\ I • ■ ■ I {-*-ml r -%m2 r • • ■ / -^mn}} 

options (List): Options to cdd (e. g. adjacency, nondegen- 

erate, minindex,...) - see documentation to cdd 

( http://www.ifor.math. ethz.ch^fukuda/cddman/cddman. html ) 

particles (Integer): Number of particles 

measurements (Integer): Number of possible measurements to each particle (equivalent 
to number of detection angles) 

Description: Creates a file with ".ext "-extension that contains the data of the given 
configuration to use in cdd. Eighter a list of vertices of the considered correlation 
polytop or the number of particles used and the possible measurements to each can be 
handed over. In the latter case the list of vertices is generated automatically. 

Example: 

in[]= ToCddExtFile ["test",2,3] 

Out [ ] = test . ext 

writes the file "test. ext" to the current working directory, containing the vertices 
of the 2-particles 3 -measurements configuration. 

3.3 TruthTable 

TruthTable [ particles, measurements, opt: options] 
particles (Integer): Number of particles 

measurements (Integer): Number of possible measurements to each particle (equivalent 
to number of detection angles) 

opt.options: The only possible option is IncludeVars — > True/False. If 
IncludeVars — * True is defined, the function includes a list of variables 
belonging to the given configuration as titles of the columns and output will be 
made in MatrixForm, otherwise a list containing all vertices is returned. De- 
fault is IncludeVars — > True. 
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Description: Creates a truth table of the given configuration, containing all vertices 
of the corresponding correlation polytopes. For generating this table all possible single 
events are rated eighter or 1 (i. e. true or false) and the joint events are evaluated 
using the logical AND operation. 

Example: 

in[]= TruthTable [2,2] 
Out[] = 
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3.4 RunCdd 

RunCdd [ file] 

file (String): File handed over to cdd as command parameter (automatically extended 
with "ext. "-suffix. 

Description: RunCdd [...] executes cdd (with file as parameter) and returns the 
corresponding H-representation, which can be used in various other functions like 
Plot Inequalities [. . . ] or Get Viol Inequalities [. . . ] . 
Using this function one has to pay attention to the potentially long runtime in calculat- 
ing the faces (i. e. the inequalities) of the correlation polytope. It can be more beneficial 
to use ToCddExtFile [ file,...] to create a ".exf"-file, followed by executing cdd 
outside of Mathematica (eventually on a faster computer) to convert the date to H- 
representation stored in an ".ine" '-file. Afterwards one can read in this file utilizing 
ReadlnHRep [ file] 
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Example: 



in[]= RunCdd [ "test"] 

out [ ] = {{H- represent at ion}, {begin}, {684, 16, real}, {2, 0,-2, . . . .}, . . . ., {end}}, 
whereas 2-particles 3-measurement configuration is taken into consideration 
here. 

3.5 ShowVRep 

ShowVRep [ particles, measurements] 
particles (Integer): Number of particles 

measurements (Integer): Number of possible measurements to each particle (equivalent 
to number of detection angles) 

Description: Shows the V-representation of a given configuration. 

Example: 

in[]= ShowVRep [ 2,2] 
Out[]= {v-representation 
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3.6 ConvToHRep 



ConvToHRep [ particles, measurements, opt: file, opt: options] or 
ConvToHRep [ vertices, opt: file, opt: options ] 

particles (Integer): Number of particles 

measurements (Integer): Number of possible measurements to each particle (equivalent 
to number of detection angles) 

vertices (List): List of m vertices in n dimensions of the form 

{{■^11/-^12/ ■ ■ • / -^lnj / {-^21 1 • • • r %2nJ r • ■ • / {-^ml r %m2 r • • • /-^/nn}}" 

file (String): Filename that is used for the conversion from a ".ext "-file to a ".ine"-Sle 
which is equivalent to a conversion from V-representation to H-representation. 
Default is "tmp " . 

options (List): Options to cdd (e.g. adjacency, nonde- 

generate, minindex,...) - see documentation to 

cdd( http://www.ifor.math. ethz.ch^fukuda/cddman/cddm an. html) 



Description: This function converts a given configuration (n particles, m measure- 
ments) or a given list of vertices from V-representation into a H-representation. As 
above in (3.4) the potentially long calculation time has to be taken into consideration, 
depending on the complexity of the problem. 



Example: 

in[]= ConvToHRep [2,3 "2.3"] 

out[]= {{H-representation}, {begin}, {684, 16, real}, {2, 
0, . . .}. . .,{end}}, 

wheras in this case the files "23. ext " (created by Mathematica containing the 
data for cdd) and "23.ine" (created by cdd as result of the calculation) are 
generated in the current working directory. 

3.7 ReadlnHRep 

ReadlnHRep [ file] 

file (String): ".ine " - file containing the H-representation which is to be read in. 



Description: Reads the H-representation from a given ".ine "-file for fur- 
ther use in various functions like GetViolInequalities [...] or 
Plotlnequalities [...] . 
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Example: 

m[]= ReadlnHRep [ "2.3"] 

out [ ] = {{H-representation},{begin},{684, 16,real},{2, 0,...}...,{end}} 

3.8 GetlnequFromHRep 

GetlnequFromHRep [hrep] 

hrep (List): H-representation yielded for example from the function 

ReadlnHRep [. . . ] or ConvToHRep [. . . ] . 

Description: Returns the inequalities from a given H-representation as a list. To 
make the list more readable one can apply InequToRead [ . . . ] on it. 

Example: 

in[]= GetlnequFromHRep [ConvToHRep [2,1] ] 

Out[]= {{al - albl + bl, l}, {-al + albl, 0}, {albl - bl, 
0}, {-albl, 0}} 

3.9 InequToRead 

InequToRead [ inequalities ] 

inequalites (List): List of inequalities yielded from GetlnequFromHRep [. . . ] 

Description: Makes the list of inequalities yielded from 

GetlnequFromHRep [...] more readable. 

Example: 

in[]= InequToRead [GetlnequFromHRep [ConvToHRep [2,1] ] ] 

Out[]= a\—a\b\+b\<\ 
-al+albl<0 
albl-bl<0 
-albl <0 

3.10 GetViolInequalities 

GetViolInequalities [ hrep, angles, functions, inequ-nr, violation] or 
GetViolInequalities [ file, angles, functions, inequ-nr, violation, opt: options] 

hrep (List): H-representation yielded for example from the function 

ReadlnHRep [. . . ] or ConvToHRep [. . . ] . 
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file (String): File containing the demanded H-representation in cdd format. 

angles (List): List of measurement angles for each particle, whereas the dimension of 
the list must represent the configuration. If one chooses the configuration "2 
particles - 2 measurements" the list must have the dimension [2,2], in this 
case the particles a and b are measured each along two axis given by the angles 
ai, «2, b\ and bi, so this parameter has the form {{a\ , a{\ , {b\ , b{\\ . 

functions (Symbol): Functions to calculate the quantum mechanical probability of the 
events. Considering for example two spin-^ particles in a singlett state, the prob- 
ability to find the particles both either in spin "up" or both in spin "down" states 
is given by P^(8,(j)) = P^(8,(j)) = | sin 2 [^^], where and (j) are the measure- 
ment angles of the particles. 

In defining these functions one has to notice, that for all possible events (sin- 
gle events, two-particle-events,. . . ) an apropriate function definition has to exist, 
each taking a list as parameter (e.g. Prob [{x_, y_} ] =jSin [ (x-y) / 2] 2 ). 

inequ-nr (List): Range of rows in H-representation used for checking violated inequal- 
ities. Specifying this can be useful, if many inequalities have to be evaluated. 
The form of the parameter is {min, max} oder All. 

violation (Real): Only inequalities are printed out, that are violated more than this pa- 
rameter. Default value is 0. 

opt.options: Options can be Printout — > True/False, which specifies, if the 
inequalities are printed out during evaluation or not. 

Description: GetViolInequalities [... ] calculates the discrepancy of in- 
equalities using the given probability functions and therefore the quantum mechanical 
violation for a distinct adjustment (i. e. special angles) of the detectors and returns all 
violated inequalities. 

Example: 

m[]= GetViolInequalities [ConvToHRep [2,2] ,{{-|,0},{0,|}},Prob] 

Out[]= (-albl+alb2+a2-a2bl-a2b2+bl<l \) 

The probability for the single event of measuring one particle in spin "up" or 
spin "down" at any angle is given by Prob [ {x_} ] =j and the probability of the 
joint event has been calculated by Prob [{x_, y_}] =jSin [ (x - y ) /2 ] 2 . 

3.11 Plotlnequalities 

Plotlnequalities [ hrep, range, angles, functions, opt: options] or 
Plotlnequalities [file, range , angles , functions , inequ-nr, violation , opt: options ] 

hrep (List): H-representation yielded for example from the function 

ReadlnHRep [. . . ] or ConvToHRep [. . . ] . 
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fiJe (String): File containing the demanded H-representation in cdd format. 

range (List): Parameter specifying the range for the variable x, which is plotted on 
the horizontal axis. It has the form {x,x m! „ , x mc!JC }(see Mathematica function 
Plot [...]) 

angles (List): List of measurement angles for each particle, whereas the dimension of 
the list must represent the configuration. If one chooses the configuration "2 
particles - 2 measurements" the list must have the dimension [2,2], in this 
case the particles a and b are measured each along two axis given by the angles 
fli, a2, b\ and bi, so this parameter has the form {{a\ , a{\ , {b\ , b{\\ . 

functions (Symbol): Functions to calculate the quantum mechanical probability of the 
events. Considering for example two spin- ^ particles in a singlett state, the prob- 
ability to find the particles both either in spin "up" or both in spin "down" states 
is given by P^(9,(j)) = P^(8,(j)) = | sin 2 [^^], where 8 and (j) are the measure- 
ment angles of the particles. 

In defining these functions one has to notice, that for all possible events (sin- 
gle events, two-particle-events,. . . ) an apropriate function definition has to exist, 
each taking a list as parameter (e.g. Prob [{x_, y_}] = jSin[ (x-y) / 2] 2 ). 

inequ-nr (List): Range of rows in H-representation used for checking violated inequal- 
ities. Specifying this can be useful, if many inequalities have to be evaluated. 
The form of the parameter is {min, max} oder All. 

violation (Real): Only inequalities are printed out, that are violated more than this pa- 
rameter. Default value is 0. 

opt: options: Options for the Mathematica- function Plot [...] can be handed over. 

Description: This function yields a plot of (violated) inequalities, whereas the func- 
tion plotted is f(x) — p{x) — b derived from the inequalites of the form p(x)<b (p(x) is 
a linear combination of functions to calculate the probabilities of (joint) events, depen- 
dent on one variable). Consequently the degree of violation is represented as a positive 
value of f(x). 

Take for example the case "2 particles and 2 measurement directions", where the 
inequality 

-Pa\b\ + Pa\b2 + Pal - Pa2b\ - Palbl + Pbl < ^ 

appears. The probability for the single event of measuring one particle in spin "up" or 
spin "down" at any angle is given by 

Pal (x) = Pal{x) = Pbl (x) = Pbl{x) = ^ 

and the probability of the joint event can be calculated by 

1 , 

Pa\b\{x,y) =p a 2bi(x,y) =Pam(x,y) =Pa2b2(x,y) = - sin[(x - y)/2] . 
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If we define the measurement angles by 



a\ = —-+x ci2 = b\=0 bi = 2*% 

the inequality can be written as 

1 . 1 n . sinx 1.1, 7t . , 
l + -sm-( — -—x) -sin-(-- +x) < 1 

2 2 V 3 2 2 2 V 3 

The left side is dependent on x (p(x) = 1 + \ sin - x) - ^ - \ sin ±(- f 
and = 1. The function to be plotted is /(x) = p(x) — b, thus 

N , 1 . 1 , 7t . sinx 1 . 1 7t , , 
fix) = 1 + - sin-( — - —x) r-sin-(— - +x) — 1 

2 2 V 3 7 2 2 2 V 3 ; 

Example: 

in[]= Plotlnequalities [ConvToHRep [2,2] ,{x,0,7t},{{- 
f,+x,0},{0,2x}},Prob] 

Out [ ] = 




The functions to calculate the probabilities have been 
defined as Prob[{x_}] :=| for a single event and 
Prob [{x_, y_}] :=5Sin[ (x-y) /2] 2 for two-particle events. 

3.12 ContPlotlnequalities 

ContPlotlnequalities [ hrep,rangex,rangey, angles, func, ineq- 
nr, violation, opt: options ] or 

ContPlotlnequalities [ file, rangex, rangey, angles, func, ineq- 
nr, violation, opt: options ] 

hrep (List): H-representation yielded for example from the function 

ReadlnHRep [. . . ] or ConvToHRep [. . . ] . 
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rangex (List): Parameter specifying the range for the variable x, which is plotted on 
the horizontal axis. It has the form {x,x mi „ , x mcK }(see Mathematica function 
Plot [...]) 

rangex (List): Parameter specifying the range for the variable y, which is plotted on 
the vertical axis. It has the form {y ,y m(n ,ymax}(see Mathematica function 
Plot [...] ) 

angles (List): List of measurement angles for each particle, whereas the dimension of 
the list must represent the configuration. If one chooses the configuration "2 
particles - 2 measurements" the list must have the dimension [2,2], in this 
case the particles a and b are measured each along two axis given by the angles 
a\, «2, b\ and b2, so this parameter has the form {{a\ , a{\ , \b\ , bi}} ■ 

functions (Symbol): Functions to calculate the quantum mechanical probability of the 
events. Considering for example two spin-^ particles in a singlett state, the prob- 
ability to find the particles both either in spin "up" or both in spin "down" states 
is given by P^(6,(j)) = P^(9,(j)) = \ sin 2 ^ 9 ^ ], where 8 and § are the measure- 
ment angles of the particles. 

In defining these functions one has to notice, that for all possible events (sin- 
gle events, two-particle-events,. . . ) an apropriate function definition has to exist, 
each taking a list as parameter (e.g. Prob [{x_, y_}] = jSin[ (x-y) / 2] 2 ). 

inequ-nr (List): Range of rows in H-representation used for checking violated inequal- 
ities. Specifying this can be useful, if many inequalities have to be evaluated. 
The form of the parameter is {min, max} oder All. 

violation (Real): Only inequalities are printed out, that are violated more than this pa- 
rameter. Default value is 0. 

opt: options: Options for the Mathematica-function Plot [...] can be handed over. 

Description: Like the function Plotlnequalities [. . . ] 

ContPlotlnequalities [...] yields a graphical representation of the vi- 
olation of Boole-Bell type inequalities, but in this case the functions are dependant on 
two variables: f(x,y) — p(x,y) — b derived from p(x,y) < b, where p(x,y) is a linear 
combination of functions to calculate the probability for single or joint events. Like in 
the description of the Plotlnequalities [... ] -function we take the "2 particles 
- 2 measurement directions", the only difference is the selection of the measurement 
angles: 

a\=x 02 = ^1=0 b2=y 
Thus the function f(x,y) is now given by 

fl \ , 1 • x 1 . x-y 1 . y 

and can be plotted as contour plots. A higher level of violation is represented by a 
darker contour layer. 
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Example: 

in[]= ContPlotlnequalities [ConvToHRep [2,2] ,{x,0,7t},{y,0,7t},{{x,0},{0,y}},Prob,All] 

Out [ ] = 

albl - alb2 + a2 - a2bl - a2b2 + b2 <= 1 al - albl - alb2 - a2bl + a2b2 + bl <= 1 



0.5 1 1.5 2 2.5 3 



0.5 1 1.5 2 2.5 3 



The functions to calculate the probabilities have been 
defined as Prob[{x_}] :=| for a single event and 
Prob [{x_, y_}] :=:jSin[ (x-y) /2 ] 2 for two-particle events. 

3.13 Cdd 

Cdd [ values, opt: command] 

values (List): Data handed over to cdd in an input-file. Each element of the list must 
be a string and represents a row in the output-file. 

opt: command (String): Command to be executed on the generated output-file (de- 
fault is "cdd" - the default value can be changed utilizing the function 
SetCddCmd [ newcommand.String] ). Here one can specify for example "cddf" 
or "cddr". 



Description: Simple Interface to run cdd in Mathematicci. The current version can- 
not distinguish, whether cdd has produced correct output or not, so one has to pay 
attention while using this function. 

Example: 

in [ ] = Cdd [ {"H-Representation","begin","6 4 real","2 -1 0","2 -1 0","-l 1 O'Y'-l 1 
0","-l 1","4 -1 -1 O'Y'end"}] 

out [ ] = {{"*", "cdd:", "Double", "Description", "Method", "C-Code: Version", 

"0.61b", "(November", 29, "1997)"}, {"*", "Copyright", "(C)", 1996, "Komei", 
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"Fukuda,", ''fukuda@ifor.math.ethz.ch''}, {"*Input", "File:tmp.ine", "(", 6, "x", 
"4)"}, {"*HyperplaneOrder:", "LexMin"}, {"*Degeneracy", "preknowledge", 
"for", "computation:", "None", "(possible", "degeneracy)"}, {"*Vertex/Ray", 
"enumeration", "is", "chosen."}, {"*Computation", "completed", "at", 
"Iteration", 6.}, {"*Computation", "starts", "at", "Thu", "Mar", 22, "18:48:36", 
2001}, {"*", "terminates", "at", "Thu", "Mar", 22, "18:48:36", 2001}, 
{"*Total", "processor", "time", "=", 0, "seconds"}, {"*", "=", 0, "hour", 0, 
"min", 0, "sec"}, {"*FINAL", "RESULT:"}, {"*Number", "of", "Vertices", 
"=", 4, "Rays", "=", 1}, {"V-representation"}, {"begin"}, {5, 4, "real"}, {1,2, 
1, 1}, {1, 1, 1, 1}, {1, 1, 2, 1}, {1, 2, 2, 1}, {0, 0, 0, 1}, {"end"}, {"hull"}} 

4 Examples 

These following two examples are originally solved in a paper by Pitowsky and Svozil 
Jl5|], The associated Mathematica - notebooks are "3_2.nb" (three particles - 2 mea- 
surement directions) and "2_3.nb" (two particles and three measurement directions) 

4.1 Three particles and two measurement directions 

In this configuration three particles (a, b and c) are measured in detectors which can be 
switched between two angles each. Consequently there are six different propositions 
for single particle events: a\, a2, b\, £>2, c\, C2, supposing that a\ is the detection (i. e. 
the click in a counter) of the particle a in the detector set along the axis specified by 
the first angle for particle a, bi the detection of particle b at the second angle for this 
particle, and so on . . . (cf. Figure |J). If we also take two and three particle events into 
account (for example the event fl2 c i means a click in the counter for particle a at the 
second angle AND a click in the counter for particle c at the first angle), there are 26 
different events: 

a\, a2, b\, b2, c\, C2, a\b\, a\b2, a2b\, 02^2, &\C\, a\C2, a2Ci, 0202, b\c\, b\C2, b2C\, 
b2, C2,a\b\c\, a\b\C2, a\b2C\, ai&2C2, 02^1^1, 02^1^2, 02^1, Q2&2C2 

4.2 Violations of inequalities 

The truth table for this configuration can be obtained utilizing the function 
TruthTable [ 3,2] (], executing ConvToHRep [3,2] yielded the appropriate H- 
representation, but this would last quite long, due to the complexity of the correla- 
tion polytope for this setting (there are 53856 hyper-planes limiting the polytope). 
Because of this fact trying to read in the H-representation created by cdd (using 
ReadlnHRep [. . . ] ) could also result in memory resource problems by Mathemat- 
ica. 

To avoid this symptoms it is suggested to export the list of vertices and apply cdd out- 
side of Mathematica to the file containing the list of vertices (V-representation). This 
can be done by invoking ToCddExtFile [ "3.2", 3,2] , which creates a file "3 _2.ext". 

2 due to lack of space not listed here, but it can be found in the Mathematica- notebook 3J2.nb 
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Figure 3: Setting for "2 particles - 3 angles" 



This file can be handed over to cdd as parameter to get the file "3_2.ine" comprising 

the H-representation of the correlation polytope (Command: "cdd 3J2.ext"). 

Now the search for violated inequalities can begin using the function 

GetViolInequalities [...] : 

May be accepted that the functions to calculate the quantum probabilities of the (joint) 
events (Prob) have been defined by Prob[{x_}] :=|, Prob [{x_, y_}] :=j and 
Prob [{x_, y_, z_}] :=g * (1 — Sin[x + y + z]), where x, y and z are the angles used 
for detection of each particle, 

GetViolInequalities ["3_2.ine",{{0,§ },{0,§ },{0,§ }},Prob,All,0.4] 

yields: 

{{-3 al+2 albl+alblcl-4 alblc2+3 alb2-3 alb2cl- 

alb2c2+alcl+3 alc2+2 a2bl-2 a2blcl-a2blc2- 
2 a2b2+a2b2cl+3 a2b2c2+a2cl-a2c2-2 bl+ 
blcl+2 blc2+b2cl-2 b2c2-cl < 0,0.5}, 

{-2 al+2 albl+alblcl-4 alblc2+2 alb2-2 alb2cl- 

alb2c2+alcl+2 alc2+3 a2bl-3 a2blcl-a2blc2- 
2 a2b2+a2b2cl+3 a2b2c2+a2cl-2 a2c2-3 bl+ 
blcl+3 blc2+b2cl-b2c2-cl < 0,0.5}, 

{-2 al+albl+alblcl-4 alblc2+2 alb2-3 alb2cl- 
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alb2c2+2 alcl+2 alc2+2 a2bl-3 a2blcl- 

a2blc2-a2b2+a2b2cl+2 a2b2c2+a2cl-a2c2- 

2 bl+2 blcl+2 blc2+b2cl-b2c2-2 cl < 0,0.5}, 
{-2 al+2 albl+alblcl-4 alblc2+2 alb2-3 alb2cl- 

alb2c2+3 alc2+2 a2bl-3 a2blcl-a2blc2- 

2 a2b2+a2b2cl+2 a2b2c2+a2cl-2 bl+2 blcl+ 

2 blc2+b2cl-b2c2-cl-c2 < 0,0.5}, 
{-al+albl+alblcl-4 alblc2+alb2-3 alb2cl- 

alb2c2+alcl+2 alc2+a2+a2bl-2 a2blcl- 

a2blc2-2 a2b2+a2b2cl+3 a2b2c2-a2c2-bl+ 

blcl+2 blc2+b2+b2cl-2b2c2-cl < 1,1.5}, 
{-al+2 albl-alblcl-3 alblc2+alb2-4 alb2cl+ 

alb2c2+3 alcl+a2+a2bl-2 a2blcl-a2blc2- 

2 a2b2+a2b2cl+3 a2b2c2-a2c2-2 bl+2 blcl+ 

2 blc2+b2+2 b2cl-3 b2c2-2 cl+c2 < 1,1.5}, 
{-2 al+2 albl-alblcl-3 alblc2+alb2-2 alb2cl- 

alb2c2+2 alcl+2 alc2+a2+a2bl-4 a2blcl+ 
a2blc2-2 a2b2+a2b2cl+3 a2b2c2+2 a2cl-3 a2c2- 
bl+3 blcl+b2-b2c2-2 cl+c2 < 1,1.5} 

{...}...} 

4.3 Graphical representation 

Using Plotlnequalities [. . . ] a graph can be created showing the violation of 
inequalities dependent on one variable. Defining the probability functions as above, 
executing 

Plotlnequalities [ "3J2.me",{x,0,n},{{0,x},{0,x},{0,x},Prob,{10000,2m00},0A] 

yields the following plot (cf. Figure |]), whereas the corresponding H-representation 
has to be stored in the file "3_2.ine", {10000,20000} indicates the range of row num- 
bers taken for calculating the graph and . 4 is the minimal degree of violation to 
include the inequality in the graph: 

To display inequalities dependent on two variables the function 
ContPlotlnequalities [...] is provided. This function shows the vio- 
lation as a contour plot, a more violated set of detection angles results in a darker 
region in the plot. 

ContPlotlnequalities [ '3.2.me'\{x,0,n},{x,0,n},{{0,x},{0,y},{x,y}},Prob,{10000,20000},0A] 

returns ContourGraphics-objects, which can be displayed for example by execut- 
ing 
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Figure 4: PlotInequalities["3^.ine",{x,0,Ti},{{0,x},{0,x},{0,x},Prob,{10000,20000},0.4] 



Show [GraphicsArray [Partition [ coat, 3[[{1, 2}, All] ] ] ] ] // 

TableForm 

(see figure ^). 

4.4 Two particles and three measurement directions 

In the case of two particles (a and b) with three properties (whereas the properties are 
three different angles of the detectors for each particle denoted by a\, ai, 03, b\, bi, b?, 
- see figure |6|) 15 different events can be found: 

{a\, a2, 03, bi, b2, b% c\, ci, C3, a\b\, a\b2, aibj,, a2b\, aibi, dibj, aj,b\, a$b2, 03^3} 



4.5 Violations of inequalities 

Using TruthTable [2,3 ] all vertices of the corresponding correlation polytope can 
be found - we get a dimension of 15 and 64 vertices as result (table ||). 
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Table 3: Truth table for 6 propositions 



Executing hrep=ConvToHRep [2,3] a H-representation of the polytope de- 
scribed by the truth table above can be created. This results in 684 hyper-planes re- 
spectively 684 inequalities from the 64 vertices limiting the polytope: 

{{H-representation}, {begin}, {684 , 16, real}, 

{2, 0, -2, 1, -1, 0, -1, 1, -1, 0, 1, 1, 1, -1, -1, 1}, 

{2, -2, 0, 1, -1, 0, -1, 1, 1, 1, 1, -1, 0, -1, -1, 1}, 

{•••}, 

{1, -1, 0, 0, -1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0}, 
{end}, {Konf iguration, 2, 3}} 
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-al + albl - 2 alblcl - 2 alblc2 + 3 alb2 - 
2 alb2cl - 2 alb2c2 + alcl + alc2 + a2bl + 
a2blcl - 3 a2blc2 - a2b2 + 3 a2b2cl - a2b2c2 - 
a2cl + a2c2 - bl + 2 blc2 - 2 b2 + 2 b2c2 - 
<= 

3 




0.5 



0.5 1 1.5 2 2.5 3 



. - 3 alblcl + alblc2 + 3 alb2 - alb2cl - 
4 alb2c2 + 2 alcl + alc2 - a2 + a2bl + 
a2blcl - 3 a2blc2 + 2 a2b2cl - a2b2c2 - a2cl + 
2 a2c2 + blcl - 2 b2 + 3 b2c2 - cl - c2 <= 
3 



2.5 




0.5 

o [ 

0.5 1 1.5 2 2.5 3 



-2 al + 3 albl - 2 alblcl - 3 alblc2 - alb2 + 

2 alb2cl - alb2c2 + alcl + 2 alc2 + a2 - 
a2bl + 2 a2blcl + a2blc2 - a2b2 - 2 a2b2cl + 

3 a2b2c2 + a2cl - 2 a2c2 - blcl + b2 + b2cl - 
2 b2c2 - cl + c2 <= 1 




1.5 



0.5 

o l 

0.5 1 1.5 2 2.5 3 



-al + 2 albl - 3 alblcl - alblc2 + alb2 - 

alb2cl - 2 alb2c2 + 2 alcl + alc2 + a2bl + 
a2blcl - 4 a2blc2 - 2 a2b2 + 3 a2b2cl + 
a2b2c2 - a2cl +■ a2c2 - 2 bl + 2 blcl + 
2 blc2 - b2cl + b2c2 - cl - c2 <= 
3 




0.5 

o l 

0.5 1 1.5 2 2.5 3 



Figure 5: ContPlotInequalities["3^.ine",{x,0,Jt},{x,0,Jt},{{0,x},{0,y},{x,y}},Prob,{10000,20000},0.4] 
Show[GraphicsArray[Partition[cont,3][[{l,2},All]] // TableForm 
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All inequalities can be displayed by 

Get InequFromHRep [hrep] // InequToRead 
The result of this operation is 

-albl + alb2 + 2 a2 - a2bl - a2b2 - a2b3 - a3 + a3bl + a3b2 - a3b3 + bl + b3 < 2 
2al - albl - alb2 - alb3 - a2bl + a2b2 - a3 + a3bl + a3b2 - a3b3 + bl + b3 < 2 
al - albl - alb3 - a2bl + a2b2 + a2b3 - a3 + a3bl + a3b2 - a3b3 + bl - b2 < 1 



bl < 1 



Using 



GetViolInequalities [hrep,{{0,f ,^},{{0,^,^}},Prob,All] // 



TableForm 



0, «2 = ^2 = x an d fl 3 = ^3 = "t taking the functions Prob[{x_}] :=j and 
Prob [{x_ , y_}] : =Sin[^f-]/2 to calculate the quantum probabilities, which is equiv- 



all inequalities can be displayed that are violated at the specific angles a\ 

M 

3 

alent to the probability to find two spin-j particles in a singlet state (|\|/) = ^jd Tl 
) — I IT))) both either in spin "up" or both in spin "down". 

-al - albl + alb2 + alb3 - a2 + a2bl + a2b3 + a3bl + a3b2 - a3b3 - bl - b2 < 

-al - albl + alb2 + alb3 + a2bl - a2b2 + a2b3 - a3 + a3bl + a3b2 - bl - b3 < 

-al + alb2 + alb3 - a2 + a2bl - a2b2 + a2b3 + a3bl + a3b2 - a3b3 - bl - b2 < 

-albl + alb2 + alb3 - a2 + a2bl - a2b2 + a2b3 - a3 + a3bl + a3b2 - b2 - b3 < 

-al + alb2 + alb3 + a2bl - a2b2 + a2b3 - a3 + a3bl + a3b2 - a3b3 - bl - b3 < 

-albl + alb2 + alb3 - a2 + a2bl + a2b3 - a3 + a3bl + a3b2 - a3b3 - b2 - b3 < 

-al + alb2 + alb3 - a2b2 + a2b3 - b3 < 

-albl + alb3 - a2 + a2bl + a2b3 - b3 < 

-albl + alb2 - a3 + a3bl + a3b2 - b2 < 

-al + alb2 + alb3 + a3b2 - a3b3 - b2 < 

a2bl - a2b2 - a3 + a3bl + a3b2 - bl < 

-a2 + a2bl + a2b3 + a3bl - a3b3 - bl < 



4.6 Graphical representation 

Like in the configuration "three particles and two angles" described above a graphical 
representation can be generated either dependent on one or dependent on two variables. 
On the one hand side by using 

Plot Inequalities [hrep,{x,0,^},{{0,x,|^},{{0,x,^}},Prob] 

we get a plot of all violated inequalities (cf. Figure 0), 
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Figure 7: Plotlnequalities[hrep,{x,0,7t},{{0,x,f },{{0,x,f }},Prob] 



on the other hand side contour plots of all inequalities violated for example more than 
0.2 can be generated by executing 

cont=ContPlot Inequalities [hrep,{x,0,7c},{y,0,lt},{{0,x,y},{{0,x,y}},Prob,All,0.2] 

To display the outcome of the calculation entering Show /@ cont results in Con- 
to urPlots of the following form (cf. Figure [$]): 
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al - albl + alb2 - alb3 + a2bl - a2b3 + 2 a3 - alb2 - alb3 + a2 + a2bl - a2b2 - a2b3 + 2 a3 - 
a3bl - a3b2 - a3b3 + bl + 2 b3 <= 3 a3bl - a3b2 - a3b3 + b2 + 2 b3 <= 3 



0.5 1 1.5 2 2.5 3 



0.5 1 1.5 2 2.5 3 



al - albl - alb2 + alb3 + a2 - a2bl - a2b3 + al - albl + alb2 - alb3 + a2bl - a2b2 - a2b3 + 



a3bl - a3b2 - a3b3 + bl + b2 <= 2 a3 - a3bl - a3b2 + bl + b3 <- 2 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



Figure 8 : ContPlotlnequalities [hrep, {x,0,7t}, {y,0,7t},{ {0,x,y }, {{0,x,y } },Prob,All,0.2] 
Show l@ % 
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